Genomics Of Drug Sensitivity In Cancer Research Articles

CDK4/6 inhibition sensitizes MEK inhibition by inhibiting cell cycle and proliferation in pancreatic ductal adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) is not sensitive to most chemotherapy drugs, leading to poor chemotherapy efficacy. Recently, Trametinib and Palbociclib have promising prospects in the treatment of pancreatic cancer. This article aims to explore the effects of Trametinib on pancreatic cancer and address the underlying mechanism of resistance as well as its reversal strategies. The GDSC (Genomics of Drug Sensitivity in Cancer) and CTD2 (Cancer Target Discovery and Development) were utilized to screen the potential drug candidate in PDAC cell lines. The dose-increase method combined with the high-dose shock method was applied to induce the Trametinib-resistant PANC-1 and MIA PaCa-2 cell lines. The CCK8 proliferation assay, colony formation assay, flow cytometry, and western blot were conducted to verify the inhibitory effect of Trametinib and Palbociclib. RNA-seq was performed in resistant PDAC cell lines to find the differential expression genes related to drug resistance and predict pathways leading to the reversal of Trametinib resistance. The GDSC and CTD2 database screening revealed that Trametinib demonstrates a significant inhibitory effect on PDAC. We found that Trametinib has a lower IC50 than Gemcitabine in PDAC cell lines. Both Trametinib and Gemcitabine can decrease the proliferation capacity of pancreatic cells, induce cell cycle arrest, and increase apoptosis. Simultaneously, the phosphorylation of the AKT and ERK pathways were inhibited by the treatment of Trametinib. In addition, the RNA-seq of Trametinib-induced resistance PDAC cell lines reveals that the cyclin-dependent kinase (CDK)-RB-E2F regulatory axis and G2/M DNA damage checkpoint might lead the drug resistance. Besides, the combination of Trametinib with Palbociclib could inhibit the proliferation and cell cycle of both resistant cells lines and also restore the sensitivity of drug-resistant cells to Trametinib. Last but not least, the interferon-α and interferon-γ expression were upregulated in resistance cell lines, which might lead to the reversal of drug resistance. The study shows Trametinib has a critical inhibitory effect on PDAC. Besides, the combination of Trametinib with Palbociclib can inhibit the proliferation of PDAC-resistant cells.

Identification of fatty acids synthesis and metabolism-related gene signature and prediction of prognostic model in hepatocellular carcinoma.

Fatty acids synthesis and metabolism (FASM)-driven lipid mobilization is essential for energy production during nutrient shortages. However, the molecular characteristics, physiological function and clinical prognosis value of FASM-associated gene signatures in hepatocellular carcinoma (HCC) remain elusive. The Gene Expression Omnibus database (GEO), the Cancer Genome Atlas (TCGA), and International Cancer Genome Consortium (ICGC) database were utilized to acquire transcriptome data and clinical information of HCC patients. The ConsensusClusterPlus was employed for unsupervised clustering. Subsequently, immune cell infiltration, stemness index and therapeutic response among distinct clusters were decoded. The tumor immune dysfunction and exclusion (TIDE) algorithm was utilized to anticipate the response of patients towards immunotherapy, and the genomics of drug sensitivity in cancer (GDSC) tool was employed to predict their response to antineoplastic medications. Least absolute shrinkage and selection operator (LASSO) regression analysis and protein-protein interaction (PPI) network were employed to construct prognostic model and identity hub gene. Single cell RNA sequencing (scRNA-seq) and CellChat were used to analyze cellular interactions. The hub gene of FASM effect on promoting tumor progression was confirmed through a series of functional experiments. Twenty-six FASM-related genes showed differential expression in HCC. Based on these FASM-related differential genes, two molecular subtypes were established, including Cluster1 and Cluster2 subtype. Compared with cluster2, Cluster1 subtype exhibited a worse prognosis, higher risk, higher immunosuppressive cells infiltrations, higher immune escape, higher cancer stemness and enhanced treatment-resistant. PPI network identified Acetyl-CoA carboxylase1 (ACACA) as central gene of FASM and predicted a poor prognosis. A strong interaction between cancer stem cells (CSCs) with high expression of ACACA and macrophages through CD74 molecule (CD74) and integrin subunit beta 1 (ITGB1) signaling was identified. Finally, increased ACACA expression was observed in HCC cells and patients, whereas depleted ACACA inhibited the stemness straits and drug resistance of HCC cells. This study provides a resource for understanding FASM heterogeneity in HCC. Evaluating the FASM patterns can help predict the prognosis and provide new insights into treatment response in HCC patients.

Abstract LB413: A newly developed patient-derived culture panel identifies precision therapies and pan-effective agents for head and neck squamous cell carcinoma

Abstract Head and neck squamous carcinoma (HNSCC) is an aggressive cancer with few treatment options. Its well-known genetic heterogeneity presents a major challenge for precision medicine development. We report a newly developed panel of 20 HNSCC patient-derived cultures (HNSCC-PDCs) from advanced or recurrent HNSCC, with comprehensive somatic and germline characterization for precision medicine (PM) discovery. The PDCs retain genetic heterogeneity of the original patient tumors, and xenografts of which also phenocopied HNSCC patient tumors. Importantly, somatic and germline events that were never or seldom captured by long-established HNSCC cell lines are represented by the PDC panel, including HRAS-CASP8 co-mutations, EGFR-AS1 germline mutation, exonic germline PIK3CA mutation, sonic hedgehog (SHH) pathway and FGFR pathway mutations. We hypothesize that drug screening with this newly developed HNSCC-PDC panel that recapitulates HNSCC patient genetics may identify new gene-drug sensitivity relationships for precision medicine development in advanced HNSCC. Here, data from our pilot drug screen of 38 selected FDA approved/Phase I-II clinical trial agents first reveal: 1) HNSCC-PDC with the hotspot HRASp.G13R, together with CASP8p.L183S, is sensitive to tipifarnib vs. tipifarnib-resistance in a long-established HNSCC cell line, SNU899 with the same HRAS mutation. Subsequent high-content imaging cell death assay further demonstrates a dose-dependent tipifarnib-induced cell death in HRASp.G13R-mutated PDC, but not in SNU-899; 2) TP53-EP300 co-mutated HNSCC-PDCs are most sensitive to cisplatin with ~340 nM sensitivity (mean IC50) among all PDCs tested, which is 22 times more sensitive than that reported in Genomics of Drug Sensitivity in Cancer (GDSC) HNSCC cell lines; 3) EGFR-AS1 (c.2361G>A)-germline mutated PDCs, and MAPK1-mutated PDCs, are both hypersensitive to EGFR inhibitors (EGFRi: erlotinib and mobocertinib), which independently credentialed gene-drug sensitivity relationships reported in exceptional responders in HNSCC clinical trials. Further signaling investigation in PDC models demonstrates phospho-MAPK activation as a reliable surrogate biomarker for EGFRi sensitivity in HNSCC; 4) FRS3-mutated PDCs are sensitive to FDA-approved FGFR pathway inhibitors; and 5) TP53-wildtype HNSCC-PDCs are 150-fold more sensitive to novel p53-MDM2 interaction inhibitors vs. TP53-mutated PDCs, a phenomenon that cannot be captured in long-established TP53-wildtype HNSCC cell lines. Lastly, 4 pan-PDC effective agents are identified with nM sensitivities and high drug activities (i.e. drug-induced growth inhibition) in more than half of our PDC models. These include 3 FDA-approved agents, TAK-981, Selinexor (a novel XPO1 inhibitor), and trametinib, as well as the preclinical agent JQ1. In conclusion, our pilot HNSCC-PDC drug screen identifies new gene-drug sensitivity relationships for HNSCC PM development and new pan-effective agents for advanced HNSCC. Citation Format: Yuen Keng Ng, Hui Li, Wenying Piao, Cecilia Pik Yuk Lau, Peony Hiu Yan Poon, Tin Yu Samuel Law, Jason Ying Kuen Chan, Yuxiong Su, Chin Wang Lau, Zigui Chen, Laura Ren Huey Ip, Thomas Chung-Kit Yeung, Sai Fung Yeung, Jason Hoi-Lam Ngan, Zoey Wing Yee Liu, Noah R. Wiese, Vivie Vo, Zoey Patel, Jeremia Onyekachi, Laili Afzali, Jasmine Zohal Afzali, Jack Nimitz Smith, Vivian Wai Yan Lui. A newly developed patient-derived culture panel identifies precision therapies and pan-effective agents for head and neck squamous cell carcinoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 2 (Late-Breaking, Clinical Trial, and Invited Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(7_Suppl):Abstract nr LB413.

DBDNMF: A Dual Branch Deep Neural Matrix Factorization method for drug response prediction.

Anti-cancer response of cell lines to drugs is in urgent need for individualized precision medical decision-making in the era of precision medicine. Measurements with wet-experiments is time-consuming and expensive and it is almost impossible for wide ranges of application. The design of computational models that can precisely predict the responses between drugs and cell lines could provide a credible reference for further research. Existing methods of response prediction based on matrix factorization or neural networks have revealed that both linear or nonlinear latent characteristics are applicable and effective for the precise prediction of drug responses. However, the majority of them consider only linear or nonlinear relationships for drug response prediction. Herein, we propose a Dual Branch Deep Neural Matrix Factorization (DBDNMF) method to address the above-mentioned issues. DBDNMF learns the latent representation of drugs and cell lines through flexible inputs and reconstructs the partially observed matrix through a series of hidden neural network layers. Experimental results on the datasets of Cancer Cell Line Encyclopedia (CCLE) and Genomics of Drug Sensitivity in Cancer (GDSC) show that the accuracy of drug prediction exceeds state-of-the-art drug response prediction algorithms, demonstrating its reliability and stability. The hierarchical clustering results show that drugs with similar response levels tend to target similar signaling pathway, and cell lines coming from the same tissue subtype tend to share the same pattern of response, which are consistent with previously published studies.

Identification and Functional Characterization of Essential Genes Related to Gefitinib Sensitivity in Lung Adenocarcinoma.

This study aimed to determine the molecular markers related to gefitinib sensitivity for guiding the prognosis of lung adenocarcinoma (LUAD) and providing new evidence for promoting the precise treatment of LUAD. Lung adenocarcinoma (LUAD) is a prevalent lung cancer subtype with inferior survival outcomes. However, gefitinib is the first molecular targeted drug approved by Food and Drug Administration (FDA) to treat advanced LUAD. Gefitinib sensitivity-related gene targets for LUAD are rarely studied. This study was designed to probe the potential molecular markers related to the sensitivity of gefitinib in LUAD. The gene expression profiles of LUAD cells in the Genomics of Drug Sensitivity in Cancer (GDSC) database were used for Weighted Gene Co-expression Network Analysis (WGCNA) to select the modules most related to gefitinib sensitivity. The Cancer Genome Atlas (TCGA) database was used to compare the expression of LUAD and para-cancerous tissues. Differentially expressed genes (DEGs) were then filtered and intersected with the highly linked genes in the module relevant to gefitinib sensitivity. Univariate Cox regression analysis was conducted to identify prognostically related genes to LUAD. The correlation between genes and drug IC50 was calculated by Spearman correlation analysis. Quantitative RT-PCR and immunofluorescence detection validated hub genes FAM13B and PFKP expressions. Among the 10 modules divided by WGCNA, the module with the most significant positive correlation and the most significant negative correlation with gefitinib sensitivity were found. FAM13B, PFKP, FGD3, RNASE1, MUC16, GJB5, and GJB3 were hub genes related to gefitinib sensitivity in LUAD. Significantly, the low expressed FAM13 in LUAD tissues positively correlated with immune response. At the same time, overexpressed PFKP in the LUAD cohort was related to an unfavorable prognosis, cell proliferation, and cell cycle. We also found that FAM13B and PFKP expressions were enhanced in LUAD cell lines. This study identified 7 critical genes related to gefitinib sensitivity in LUAD. Functionally, genes positively correlated with gefitinib sensitivity might regulate the progression of LUAD through the immune, cell cycle, and metabolic pathways and showed potential effects in predicting sensitivity to different drugs. These findings help offer a theoretical direction for personalized treatment of LUAD.

Conserved methylation signatures associate with the tumor immune microenvironment and immunotherapy response

BackgroundAberrant DNA methylation is a major characteristic of cancer genomes. It remains unclear which biological processes determine epigenetic reprogramming and how these processes influence the variants in the cancer methylome, which can further impact cancer phenotypes.MethodsWe performed pairwise permutations of 381,900 loci in 569 paired DNA methylation profiles of cancer tissue and matched normal tissue from The Cancer Genome Atlas (TCGA) and defined conserved differentially methylated positions (DMPs) based on the resulting null distribution. Then, we derived independent methylation signatures from 2,465 cancer-only methylation profiles from the TCGA and 241 cell line-based methylation profiles from the Genomics of Drug Sensitivity in Cancer (GDSC) cohort using nonnegative matrix factorization (NMF). We correlated DNA methylation signatures with various clinical and biological features, including age, survival, cancer stage, tumor immune microenvironment factors, and immunotherapy response. We inferred the determinant genes of these methylation signatures by integrating genomic and transcriptomic data and evaluated the impact of these signatures on cancer phenotypes in independent bulk and single-cell RNA/methylome cohorts.ResultsWe identified 7,364 differentially methylated positions (2,969 Hyper-DMPs and 4,395 Hypo-DMPs) in nine cancer types from the TCGA. We subsequently retrieved three highly conserved, independent methylation signatures (Hyper-MS1, Hypo-MS1, and Hypo-MS4) from cancer tissues and cell lines based on these Hyper and Hypo-DMPs. Our data suggested that Hypo-MS4 activity predicts poor survival and is associated with immunotherapy response and distant tumor metastasis, and Hypo-MS4 activity is related to TP53 mutation and FOXA1 binding specificity. In addition, we demonstrated a correlation between the activities of Hypo-MS4 in cancer cells and the fractions of regulatory CD4 + T cells with the expression levels of immunological genes in the tumor immune microenvironment.ConclusionsOur findings demonstrated that the methylation signatures of distinct biological processes are associated with immune activity in the cancer microenvironment and predict immunotherapy response.

Exploration of Prognostic and Treatment Markers in Hepatocellular Carcinoma via GPCR-Related Genes Analysis

BackgroundG protein-coupled receptors (GPCRs), the biggest family of signaling receptors, account for 34% of all the drug targets approved by the Food and Drug Administration (FDA). It has been gradually recognized that GPCRs are of significance for tumorigenesis, but in-depth studies are still required to explore specific mechanisms. In this study, the role of GPCRs in hepatocellular carcinoma (HCC) was elucidated, and GPCR-related genes were employed for building a risk-score model for the prognosis and treatment efficacy prediction of HCC patients. MethodsPatients’ data on HCC were sourced from the Liver Hepatocellular Carcinoma-Japan (LIRI-JP) and The Cancer Genome Atlas (TCGA) databases, while GPCR-related genes were obtained from the Molecular Signatures Database (MSigDB). Univariant and multivariant Cox regression analyses, as well as least absolute shrinkage and selection operator (LASSO) were performed with the aim of identifying differentially expressed GPCR-related genes and grouping patients. Differential expression and functional enrichment analyses were performed; protein-protein interaction (PPI) mechanisms were explored; hub genes and micro ribonucleic acid (miRNA)-target gene regulatory networks were constructed. The tumor immune dysfunction and exclusion (TIDE) algorithm was utilized to evaluate immune infiltration levels and genetic variations. Sensitivity to immunotherapy and common antitumor drugs was predicted via the database Genomics of Drug Sensitivity in Cancer (GDSC). ResultsA GPCR-related risk score containing eight GPCR-related genes (atypical chemokine receptor 3 (ACKR3), C-C chemokine receptor type 3 (CCR3), CCR7, frizzled homolog 5 (FZD5), metabotropic glutamate receptor 8 (GRM8), hydroxycarboxylic acid receptor 1 (HCAR1), 5-hydroxytryptamine receptor 5A (HTR5A) and nucleotide-binding oligomerization domain-like receptor family pyrin domain containing 6 (NLRP6)) was set up. In addition, patients were classified into groups with high and low risks. Patients in the high-risk group exhibited a worse prognosis but demonstrated a more favorable immunotherapy response rate compared with those in the low-risk group. Distinct sensitivity to chemotherapeutic drugs was observed. A clinical prediction model on the basis of GPCR-related risk scores was constructed. Areas under the curves (AUC) corresponding to one-, three- and five-year survival were 0.731, 0.765 and 0.731, respectively. ConclusionsIn this study, an efficient HCC prognostic prediction model was constructed by only GPCR-related genes, which are all potential targets for HCC treatment.

Exploring the metastasis-related biomarker and carcinogenic mechanism in liver cancer based on single cell technology

BackgroundHepatocellular carcinoma (HCC) is a fatal primary malignancy characterized by high invasion and migration. We aimed to explore the underlying metastasis-related mechanism supporting the development of HCC. MethodsThe dataset of single cell RNA-seq (GSE149614) were collected for cell clustering by using the Seurat R package, the FindAllMarkers function was used to find the highly expression and defined the cell cluster. The WebGestaltR package was used for the GO and KEGG function analysis of shared genes, the Gene Set Enrichment Analysis (GSVA) was performed by clusterProfiler R package, the hTFtarget database was used to identify the crucial transcription factors (TFs), the Genomics of Drug Sensitivity in Cancer (GDSC) database was used for the drug sensitivity analysis. Finally, the overexpression and trans-well assay was used for gene function analysis. ResultsWe obtained 9 cell clusters from the scRNA-seq data, including the nature killer (NK)/T cells, Myeloid cells, Hepatocytes, Epithelial cells, Endothelial cells, Plasma B cells, Smooth muscle cells, B cells, Liver bud hepatic cells. Further cell ecological analysis indicated that the Hepatocytes and Endothelial cell cluster were closely related to the cancer metastasis. Subsequently, the NDUFA4L2-Hepatocyte, GTSE1-Hepatocyte, ENTPD1-Endothelial and NDUFA4L2-Endothelial were defined as metastasis-supporting cell clusters, in which the NDUFA4L2-Hepatocyte cells was closely related to angiogenesis, while the NDUFA4L2-Endothelial was related with the inflammatory response and complement response. The overexpression and trans-well assay displayed that NDUFA4L2 exhibited clearly metastasis-promoting role in HCC progression. ConclusionWe identified and defined 4 metastasis-supporting cell clusters by using the single cell technology, the specify shared gene was observed and played crucial role in promoting cancer progression, our findings were expected to provide new insight in control cancer metastasis.

Identification of biological markers in cancer disease using explainable artificial intelligence

AbstractThe research aims to improve the prediction of drug sensitivity on cancer cell lines using gene expression data and molecular fingerprints of drugs. The proposed study uses a deep learning model, BioMarkerX, trained on the Cancer Cell Line Encyclopedia (CCLE) and Genomics of Drug Sensitivity in Cancer (GDSC) datasets utilizing Particle Swarm Optimization technique to select specific genes as features. The model achieves high prediction accuracy with a Root Mean Square Error (RMSE) of 0.40 ± 0.02 and R2 of 0.83 ± 0.03 on the CCLE dataset, and an RMSE of 0.36 ± 0.05 and R2 of 0.83 ± 0.03 on the GDSC dataset. The approach also used an explainable artificial intelligence model to discover biological markers linked to cancer development. This can provide insights into targeted therapies for improving cancer treatment outcomes. Overall, the study presents an effective approach for identifying important biological markers relevant to cancer disease, aiding in the development of more efficient anticancer medications.

A hypoxia-related genes prognostic risk model, and mechanisms of hypoxia contributing to poor prognosis through immune microenvironment and drug resistance in acute myeloid leukemia.

Objective: Acute myeloid leukemia (AML) is a malignant hematologic cancer with poor prognosis. Emerging evidence suggests a close association between AML progression and hypoxia. The purpose of this study was to establish a new risk prognostic model for AML based on hypoxia-related genes, and to explore the mechanisms by which hypoxia-related genes affect the prognosis of AML based on tumor immune microenvironment (TIME) and drug resistance. Methods: The AML patient samples obtained from Therapeutically Applicable Research to Generate Effective Treatments (TARGET) database were classified into C1 and C2 based on hypoxia-related genes, followed by analysis utilizing Gene Ontology (GO), Kyoto Encyclopaedia of Genes and Genomes (KEGG) and Gene Set Enrichment Analysis (GSEA). Through univariate and LASSO Cox regression analysis, the hypoxia-related hub genes 26S proteasome non-ATPase regulatory subunit 11 (PSMD11) and 26S proteasome non-ATPase regulatory subunit 14 (PSMD14) were identified to construct the model. AML patient samples were obtained from the TARGET and The Cancer Genome Atlas (TCGA) databases, serving as the training and the validation sets, and were stratified into high-risk and low-risk group according to the median risk score. The correlations between the model and TIME and anti-tumor drugs were analysed using CIBERSORT and Genomics of Drug Sensitivity in Cancer (GDSC) databases. The expressions of PSMD11/PSMD14 in clinical samples and AML sensitive and drug-resistant cell lines were detected by Western blot and real-time PCR. Results: The C1 group with high expression of hypoxia-related genes had lower overall survival (OS). Immune-related signaling pathways were different between C1/C2, and hypoxia was positively correlated with the activation of mammalian target of rapamycin (mTOR) signaling pathway. The model had good accuracy in both the training and the validation sets. The high-risk group exhibited lower OS and TIME activity, and was more sensitive to several anti-tumor drugs. PSMD11/PSMD14 were highly expressed in relapsed patients and AML drug-resistant cell lines. Conclusion: The established novel risk prognostic model and experiment results offer valuable insights for predicting AML prognosis and guiding drug selection. It also provides a fundamental framework for the mechanisms through which hypoxia impacts AML prognosis by modulating TIME and drug resistance.

Prediction Model for Therapeutic Responses in Ovarian Cancer Patients Using Paclitaxel-Resistant Immune-Related lncRNAs.

Ovarian cancer (OC) is the deadliest malignant tumor in women with a poor prognosis due to drug resistance and lack of prediction tools for therapeutic responses to anti- cancer drugs. The objective of this study was to launch a prediction model for therapeutic responses in OC patients. The RNA-seq technique was used to identify differentially expressed paclitaxel (PTX)- resistant lncRNAs (DE-lncRNAs). The Cancer Genome Atlas (TCGA)-OV and ImmPort database were used to obtain immune-related lncRNAs (ir-lncRNAs). Univariate, multivariate, and LASSO Cox regression analyses were performed to construct the prediction model. Kaplan- Meier plotter, Principal Component Analysis (PCA), nomogram, immune function analysis, and therapeutic response were applied with Genomics of Drug Sensitivity in Cancer (GDSC), CIBERSORT, and TCGA databases. The biological functions were evaluated in the CCLE database and OC cells. The RNA-seq defined 186 DE-lncRNAs between PTX-resistant A2780-PTX and PTXsensitive A2780 cells. Through the analysis of the TCGA-OV database, 225 ir-lncRNAs were identified. Analyzing 186 DE-lncRNAs and 225 ir-lncRNAs using univariate, multivariate, and LASSO Cox regression analyses, 9 PTX-resistant immune-related lncRNAs (DEir-lncRNAs) acted as biomarkers were discovered as potential biomarkers in the prediction model. Single-cell RNA sequencing (scRNA-seq) data of OC confirmed the relevance of DEir-lncRNAs in immune responsiveness. Patients with a low prediction score had a promising prognosis, whereas patients with a high prediction score were more prone to evade immunotherapy and chemotherapy and had poor prognosis. The novel prediction model with 9 DEir-lncRNAs is a valuable tool for predicting immunotherapeutic and chemotherapeutic responses and prognosis of patients with OC.

Increased ACSL6 Expression Predicts a Favorable Prognosis in Triple-negative Breast Cancer.

Long-chain acyl-coenzyme A synthases (ACSLs) are responsible for the catalysis of fatty acids into their corresponding fatty acyl-CoAs. The dysregulation of ACSLs has been increasingly recognized in cancer patients. However, the function of ACSL6 in triple-negative breast cancer (TNBC) is still completely unknown. In this study, immunohistochemistry was applied to detect ACSL6 protein expression using a TNBC tissue microarray. Additionally, the mRNA levels of ACSL6 in human normal tissues and pancancer tissues were analyzed using Genotype Tissue Expression (GTEx) datasets and The Cancer Genome Atlas (TCGA) database. The correlations between the levels of ACSL6 expression and clinical characteristics were analyzed. The survival analysis of ACSL6 in TNBC was carried out using the Kaplan‒Meier Plotter online tool. Associations of ACSL6 with immune infiltration analyses were conducted using the ESTIMATE, CIBERSORT, and TISIDB databases. The relationship between ACSL6 and sensitivity to drugs was analyzed from Genomics of Drug Sensitivity in Cancer (GDSC). The results indicated a significant increase in ACSL6 expression in TNBC tissues compared to adjacent normal tissues. However, high ACSL6 expression was significantly associated with favorable survival outcomes in TNBC patients. Enrichment analysis revealed that coexpressed genes of ACSL6 were significantly enriched in various immunity processes. ACSL6 was positively correlated with the infiltration of memory CD4 T cells, while a negative correlation was found between ACSL6 and M2 macrophages and resting dendritic cells. Further analysis revealed that high levels of ACSL6 correlated with increased survival outcomes in cancer patients who received immunotherapy. Altogether, the current findings highlight the potential value of ACSL6 as a diagnostic and prognostic marker in the treatment of TNBC.

Multidimensional Transcriptomics Unveils RNF34 as a Prognostic Biomarker and Potential Indicator of Chemotherapy Sensitivity in Wilms' Tumour.

Nephroblastoma, colloquially known as Wilms' tumour (WT), is the predominant malignant renal neoplasm arising in the paediatric population. Modern therapeutic approaches for WT incorporate a synergistic combination of surgical intervention, radiotherapy, and chemotherapy, which substantially ameliorate the overall patient survival rate. Despite this, the optimal sequence of chemotherapy and surgical intervention remains a matter of contention, with each strategy presenting its own strengths and weaknesses that could influence clinical decision-making. To make some headway on this clinical dilemma, we deployed a multidimensional transcriptomics integration approach by analysing bulk RNA sequencing data with 136 samples, as well as single-nucleus RNA sequencing (snRNA-seq) and paired spatial transcriptome sequencing (stRNA) data from 32 WT specimens. Our findings identified a distinct elevation of RNF34 expression within WT samples, which correlated with unfavourable prognostic outcomes. Leveraging the Genomics of Drug Sensitivity in Cancer (GDSC), we simultaneously revealed that patients with high expression of RNF34 have higher sensitivity to commonly used chemotherapy drugs for WT. Furthermore, our analysis of snRNA and stRNA data unveiled a reduced proportion of RNF34 expression in neoplastic cells after chemotherapy. Moreover, stRNA data delineated a significant association between a higher proportion of RNF34 expression in cancer cells and adverse features such as anaplastic histology and tumour recurrence. Intriguingly, we also observed a close association between elevated RNF34 expression and a characteristic exhausted tumour immune microenvironment. Collectively, our findings underscore the pivotal role of RNF34 in the prognostic prediction potential and treatment sensitivity of WT. This comprehensive analysis can potentially inform and refine clinical decision-making for WT patients and guide future studies towards the development of optimized, rational therapeutic strategies.

Extensive prediction of drug response in mutation-subtype-specific LUAD with machine learning approach.

Lung cancer is the most prevalent cancer diagnosis and the leading cause of cancer death worldwide. Therapeutic failure in lung cancer (LUAD) is heavily influenced by drug resistance. This challenge stems from the diverse cell populations within the tumor, each having unique genetic, epigenetic, and phenotypic profiles. Such variations lead to varied therapeutic responses, thereby contributing to tumor relapse and disease progression. The Genomics of Drug Sensitivity in Cancer (GDSC) database was used in this investigation to obtain the mRNA expression dataset, genomic mutation profile, and drug sensitivity information of NSCLS. Machine Learning (ML) methods, including Random Forest (RF), Artificial Neurol Network (ANN), and Support Vector Machine (SVM), were used to predict the response status of each compound based on the mRNA and mutation characteristics determined using statistical methods. The most suitable method for each drug was proposed by comparing the prediction accuracy of different ML methods, and the selected mRNA and mutation characteristics were identified as molecular features for the drug-responsive cancer subtype. Finally, the prognostic influence of molecular features on the mutational subtype of LUAD in publicly available datasets. Our analyses yielded 1,564 gene features and 45 mutational features for 46 drugs. Applying the ML approach to predict the drug response for each medication revealed an upstanding performance for SVM in predicting Afuresertib drug response (area under the curve [AUC] 0.875) using CIT, GAS2L3, STAG3L3, ATP2B4-mut, and IL15RA-mut as molecular features. Furthermore, the ANN algorithm using 9 mRNA characteristics demonstrated the highest prediction performance (AUC 0.780) in Gefitinib with CCL23-mut. This work extensively investigated the mRNA and mutation signatures associated with drug response in LUAD using a machine-learning approach and proposed a priority algorithm to predict drug response for different drugs.

Correlation between IL3 signaling pathway-related genes and immune checkpoint inhibitor efficacy in patients with renal cell carcinoma.

There is a lack of effective biomarkers that predict immunotherapy efficacy in clear cell renal cell carcinoma(KIRC). We aimed to identify biomarkers that would predict the efficacy of KIRC treatment with immune checkpoint inhibitors (ICIs). Cohort data of KIRC patients with somatic mutations, mRNA expression and survival data from The Cancer Genome Atlas (TCGA) database and immunotherapy cohort and Genomics of Drug Sensitivity in Cancer (GDSC) database were analyzed and divided into interleukin 3 (IL3) pathway-related genes high expression (IL3-High) and IL3 pathway-related genes low expression (IL3-Low) groups according to pathway expression status to assess the relationship between the IL3 pathway-related genes activation status and the prognosis of KIRC patients treated with ICIs. The data were validated by immunohistochemistry experiments, and possible mechanisms of action were explored at the level of gene mutation landscape, immune microenvironment characteristics, transcriptome and copy number variation(CNV) characteristicsRESULTS: The IL3 pathway-related genes was an independent predictor of the efficacy of ICIs in KIRC patients, and the IL3-High group had a longer overall survival (OS); KIRC patients in the IL3-High group had increased levels of chemokines, cytolysis, immune checkpoint gene expression and abundant immunity. The IL3-Low group had poor immune cell infiltration and significant downregulation of complement activation, cytophagy, B-cell activation, and humoral immune response pathways. The high group was more sensitive to targeted drugs of some signaling pathways, and its efficacy in combining these drugs with immunity has been predicted in the published literature. The IL3 pathway-related genes can be used as a predictor of the efficacy of ICIs in KIRC. The IL3 pathway-related genes may affect the therapeutic efficacy of ICIs by affecting the expression of immune-related molecules, immune cell infiltration, and the level of immune response pathways.

Transcriptomic correlates of cell cycle checkpoints with distinct prognosis, molecular characteristics, immunological regulation, and therapeutic response in colorectal adenocarcinoma.

Colorectal adenocarcinoma (COAD), accounting for the most common subtype of colorectal cancer (CRC), is a kind of malignant digestive tumor. Some cell cycle checkpoints (CCCs) have been found to contribute to CRC progression, whereas the functional roles of a lot of CCCs, especially the integrated role of checkpoint mechanism in the cell cycle, remain unclear. The Genomic Data Commons (GDC) The Cancer Genome Atlas (TCGA) COAD cohort was retrieved as the training dataset, and GSE24551 and GSE29623 were downloaded from Gene Expression Omnibus (GEO) as the validation datasets. A total of 209 CCC-related genes were derived from the Gene Ontology Consortium and were subsequently enrolled in the univariate, multivariate, and least absolute shrinkage and selection operator (LASSO) Cox regression analyses, finally defining a CCC signature. Cell proliferation and Transwell assay analyses were utilized to evaluate the functional roles of signature-related CCCs. The underlying CCC signature, molecular characteristics, immune-related features, and therapeutic response were finally estimated. The Genomics of Drug Sensitivity in Cancer (GDSC) database was employed for the evaluation of chemotherapeutic responses. The aberrant gene expression of CCCs greatly contributed to COAD development and progression. Univariate Cox regression analysis identified 27 CCC-related genes significantly affecting the overall survival (OS) of COAD patients; subsequently, LASSO analysis determined a novel CCC signature. Noticeably, CDK5RAP2, MAD1L1, NBN, RGCC, and ZNF207 were first identified to be correlated with the prognosis of COAD, and it was proven that all of them were significantly correlated with the proliferation and invasion of HCT116 and SW480 cells. In TCGA COAD cohort, CCC signature robustly stratified COAD patients into high and low CCC score groups (median OS: 57.24 months vs. unreached, p< 0.0001), simultaneously, with the good AUC values for OS prediction at 1, 2, and 3 years were 0.74, 0.78, and 0.77. Furthermore, the prognostic capacity of the CCC signature was verified in the GSE24551 and GSE29623 datasets, and the CCC signature was independent of clinical features. Moreover, a higher CCC score always indicated worse OS, regardless of clinical features, histological subtypes, or molecular subgroups. Intriguingly, functional enrichment analysis confirmed the CCC score was markedly associated withextracellular, matrix and immune (chemokine)-related signaling, cell cycle-related signaling, and metabolisms. Impressively, a higher CCC score was positively correlated with a majority of chemokines, receptors, immunostimulators, and anticancer immunity, indicating a relatively immune-promoting microenvironment. In addition, GSE173839, GSE25066, GSE41998, and GSE194040 dataset analyses of the underlying CCC signature suggested that durvalumab with olaparib and paclitaxel, taxane-anthracycline chemotherapy, neoadjuvant cyclophosphamide/doxorubicin with ixabepilone or paclitaxel, and immunotherapeutic strategies might be suitable for COAD patients with higher CCC score. Eventually, the GDSC database analysis showed that lower CCC scores were likely to be more sensitive to 5-fluorouracil, bosutinib, gemcitabine, gefitinib, methotrexate, mitomycin C, and temozolomide, while patients with higher CCC score seemed to have a higher level of sensitivity to bortezomib and elesclomol. The novel CCC signature exhibited a good ability for prognosis prediction for COAD patients, and the CCC score was found to be highly correlated with molecular features, immune-related characteristics, and therapeutic responses, which would greatly promote clinical management and precision medicine for COAD.

Exploration and validation of key genes associated with early lymph node metastasis in thyroid carcinoma using weighted gene co-expression network analysis and machine learning.

Thyroid carcinoma (THCA), the most common endocrine neoplasm, typically exhibits an indolent behavior. However, in some instances, lymph node metastasis (LNM) may occur in the early stages, with the underlying mechanisms not yet fully understood. LNM potential was defined as the tumor's capability to metastasize to lymph nodes at an early stage, even when the tumor volume is small. We performed differential expression analysis using the 'Limma' R package and conducted enrichment analyses using the Metascape tool. Co-expression networks were established using the 'WGCNA' R package, with the soft threshold power determined by the 'pickSoftThreshold' algorithm. For unsupervised clustering, we utilized the 'ConsensusCluster Plus' R package. To determine the topological features and degree centralities of each node (protein) within the Protein-Protein Interaction (PPI) network, we used the CytoNCA plugin integrated with the Cytoscape tool. Immune cell infiltration was assessed using the Immune Cell Abundance Identifier (ImmuCellAI) database. We applied the Least Absolute Shrinkage and Selection Operator (LASSO), Support Vector Machine (SVM), and Random Forest (RF) algorithms individually, with the 'glmnet,' 'e1071,' and 'randomForest' R packages, respectively. Ridge regression was performed using the 'oncoPredict' algorithm, and all the predictions were based on data from the Genomics of Drug Sensitivity in Cancer (GDSC) database. To ascertain the protein expression levels and subcellular localization of genes, we consulted the Human Protein Atlas (HPA) database. Molecular docking was carried out using the mcule 1-click Docking server online. Experimental validation of gene and protein expression levels was conducted through Real-Time Quantitative PCR (RT-qPCR) and immunohistochemistry (IHC) assays. Through WGCNA and PPI network analysis, we identified twelve hub genes as the most relevant to LNM potential from these two modules. These 12 hub genes displayed differential expression in THCA and exhibited significant correlations with the downregulation of neutrophil infiltration, as well as the upregulation of dendritic cell and macrophage infiltration, along with activation of the EMT pathway in THCA. We propose a novel molecular classification approach and provide an online web-based nomogram for evaluating the LNM potential of THCA (http://www.empowerstats.net/pmodel/?m=17617_LNM). Machine learning algorithms have identified ERBB3 as the most critical gene associated with LNM potential in THCA. ERBB3 exhibits high expression in patients with THCA who have experienced LNM or have advanced-stage disease. The differential methylation levels partially explain this differential expression of ERBB3. ROC analysis has identified ERBB3 as a diagnostic marker for THCA (AUC=0.89), THCA with high LNM potential (AUC=0.75), and lymph nodes with tumor metastasis (AUC=0.86). We have presented a comprehensive review of endocrine disruptor chemical (EDC) exposures, environmental toxins, and pharmacological agents that may potentially impact LNM potential. Molecular docking revealed a docking score of -10.1 kcal/mol for Lapatinib and ERBB3, indicating a strong binding affinity. In conclusion, our study, utilizing bioinformatics analysis techniques, identified gene modules and hub genes influencing LNM potential in THCA patients. ERBB3 was identified as a key gene with therapeutic implications. We have also developed a novel molecular classification approach and a user-friendly web-based nomogram tool for assessing LNM potential. These findings pave the way for investigations into the mechanisms underlying differences in LNM potential and provide guidance for personalized clinical treatment plans.

A prognostic signature for lung adenocarcinoma by five genes associated with chemotherapy in lung adenocarcinoma.

Lung adenocarcinoma (LUAD) is one of the most common subtypes of lung cancer. Finding prognostic biomarkers is helpful in stratifying LUAD patients with different prognosis. We explored the correlation of LUAD prognosis and genes associated with chemotherapy in LUAD and obtained data of LUAD patients from the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. Drug sensitivity data were acquired from the Genomics of Drug Sensitivity in Cancer (GDSC) database. Differential and enrichment analyses were used to screen the target genes utilizing limma and "clusterProfiler" packages. Then univariate and LASSO Cox analyses were used to select the prognosis-related genes. Survival analysis was used to estimate the overall survival (OS) of different groups. Twenty-three differentially expressed genes (DEGs) were screened between LUAD samples and healthy samples, and BTK, FGFR2, PIM2, CHEK1, and CDK1 were selected to construct a prognostic signature. The OS of patients in the high-risk group (risk score higher than 0.69) was worse than that in the low-risk group (risk score lower than 0.69). The risk score model constructed by five genes is a potential prognostic biomarker for LUAD patients.

Upregulated PARP1 confers breast cancer resistance to CDK4/6 inhibitors via YB-1 phosphorylation

BackgroundCyclic-dependent kinase (CDK) 4/6 kinases, as the critical drivers of the cell cycle, are involved in the tumor progression of various malignancies. Pharmacologic inhibitors of CDK4/6 have shown significant clinical prospects in treating hormone receptor-positive and human epidermal growth factor receptor-negative (HR + /HER2-) breast cancer (BC) patients. However, acquired resistance to CDK4/6 inhibitors (CDK4/6i), as a common issue, has developed rapidly. It is of great significance that the identification of novel therapeutic targets facilitates overcoming the CDK4/6i resistance. PARP1, an amplified gene for CDK4/6i-resistant patients, was found to be significantly upregulated during the construction of CDK4/6i-resistant strains. Whether PARP1 drives CDK4/6i resistance in breast cancer is worth further study.MethodPARP1 and p-YB-1 protein levels in breast cancer cells and tissues were quantified using Western blot (WB) analysis, immunohistochemical staining (IHC) and immunofluorescence (IF) assays. Bioinformatics analyses of Gene Expression Profiling Interactive Analysis (GEPIA), Genomics of Drug Sensitivity in Cancer (GDSC) and Cancer Cell Line Encyclopedia (CCLE) datasets were applied to explore the relationship between YB-1/PARP1 protein levels and CDK4/6i IC50. Cell Counting Kit-8 (CCK-8) and crystal violet staining assays were performed to evaluate cell proliferation rates and drug killing effects. Flow cytometry assays were conducted to assess apoptosis rates and the G1/S ratio in the cell cycle. An EdU proliferation assay was used to detect the DNA replication ratio after treatment with PARP1 and YB-1 inhibitors. A ChIP assay was performed to assess the interaction of the transcription factor YB-1 and associated DNA regions. A double fluorescein reporter gene assay was designed to assess the influence of WT/S102A/S102E YB-1 on the promoter region of PARP1. Subcutaneous implantation models were applied for in vivo tumor growth evaluations.ResultsHere, we reported that PARP1 was amplified in breast cancer cells and CDK4/6i-resistant patients, and knockdown or inhibition of PARP1 reversed drug resistance in cell experiments and animal models. In addition, upregulation of transcription factor YB-1 also occurred in CDK4/6i-resistant breast cancer, and YB-1 inhibition can regulate PARP1 expression. p-YB-1 and PARP1 were upregulated when treated with CDK4/6i based on the WB and IF results, and elevated PARP1 and p-YB-1 were almost simultaneously observed during the construction of MCF7AR-resistant strains. Inhibition of YB-1 or PAPR1 can cause decreased DNA replication, G1/S cycle arrest, and increased apoptosis. We initially confirmed that YB-1 can bind to the promoter region of PARP1 through a ChIP assay. Furthermore, we found that YB-1 phosphorylated at S102 was crucial for PARP1 transcription according to the double fluorescein reporter gene assay. The combination therapy of YB-1 inhibitors and CDK4/6i exerted a synergistic antitumor effect in vitro and in vivo. The clinical data suggested that HR + /HER2- patients with low expression of p-YB-1/PARP1 may be sensitive to CDK4/6i in breast cancer.ConclusionThese findings indicated that a ‘‘YB-1/PARP1’’ loop conferred resistance to CDK4/6 inhibitors. Furthermore, interrupting the loop can enhance tumor killing in the xenograft tumor model, which provides a promising strategy against drug resistance in breast cancer.

Integrative Multi-Omic Analysis for Prognosis Stratification in Acute Myeloid Leukemia

Acute myeloid leukemia (AML) is a heterogeneous hematopoietic malignancy with a dismal prognosis. European LeukemiaNet (ELN) is crucial for tailoring AML treatment individually. Various AML models correlated survival and clinical drug response with immune cell differentiation state by deconvoluting transcriptomics. However, few comprehensive molecular subtyping models integrate multi-omic profiles for prognostic and drug response prediction. Thus, this study aimed to integrate DNA methylation, genomics, transcriptomics and ex vivo drug sensitivity screening of AML patients and explore their roles in facilitating AML molecular stratification and redefining prognosis. Firstly，data on 90 patients with AML(non-promyelocytic leukemia) were retrieved from The Cancer Genome Atlas (TCGA) for consensus clustering by 10 advanced consensus clustering algorithms. The optimal three distinct subtypes, UAMOCS1, UAMOCS2, and UAMCOS3, were identified to achieve the best clustering effect. Multi-omics were orchestrated together for UAMOCS shaping. The UAMOCS can significantly distinguish overall survival (median 18 vs. 27 vs. 45 months, P= 0.0048) among the three instinct subtypes with the integrated contribution of transcriptomic profiles and DNA methylation and mutation data. Also, the activity of AML-specific transcription factors and chromatin remodeling makers showed a diverse pattern among the three subtypes, epigenetically contributing to UAMOCS shaping. Next, the nearest template prediction, a model-free method, was applied to further validate this predictive model. Similar subtyping under this model is generated robustly by imitating 50 subtype-specific genes and capable of identifying survival in both the “ihCAMs-AML” (Figure 1) and GSE37642. The main validation cohort was a real-world de novo AML cohort (acute promyelocytic leukemia excluded), “ihCAMs-AML”（N=98）, aged 14-70 years, with a median follow-up of 28 (1-55) months. The samples were simultaneously available for methylome (N=31, 850K Methylation Chip), DNA mutation (N=96, targeted deep sequencing), transcriptome (N=98, bulk RNA-seq), and ex vivo drug screening (N=66, 17 drugs under 100% peak plasma concentration) profiling analysis. Clinical relevance showed that traditional cytogenetic stratification was strongly correlated with our system. UAMOCS1, acting as an “immune activator”, is characterized by enriching of AML myelodysplasia-related (AML-MR) gene mutations (especially RUNX1 mutation, P = 0.00839), unstable chromosome (deletions on 7q), increased abundance of immune cells with immune escaping marker and poor risk. UAMOCS2 was defined as an intermediate immune burden and a “monocytic-like” phenotype with intermediate survival, corresponding to its inclining result toward M4/M5. UAMOCS3 was an “immune desert” group supposed to be the most favorable subset among the three groups. The ex vivo drug screening result implied that UAMOCS also provide values in distinguishing drug response. UAMOCS1 hardly benefited from commonly used chemotherapeutic drugs compared with the other two distinct subtypes. Pharmacogenomics database GDSC (the Genomics of Drug Sensitivity in Cancer) among AML cell lines was used to evaluate potential drug targets for overcoming drug resistance. Receptor tyrosine kinase signaling inhibitors (SB505124 and EphB4) might be sensitive to UAMOCS1. Taken together, UAMOCS model was generated that can classify AML patients independently and help develop precise chemotherapeutic strategy for patients with AML in the future.